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Abstract: A new Monte-Carlo generator including real radiated photons in DIS 
K* ! on polarized and unpolarized targets is presented. Analytical and numerical tests 

Q I are performed and discussed in details. 

a^ ■ 1 Introduction 

p • Experimental data on lepton-nucleon scattering contain contributions from the usual Born pro- 
A cess and from QED radiative effects which come from loop corrections and from processes with 
<D , the emission of additional real photons. Any Monte-Carlo program simulating the experimen- 
.. \ tal situation has to take into account these radiative processes by generating radiated photons 
because radiative corrections can be very large. 

Depending on the four-momentum transfer squared {Q^) and the energy transfer (u) there 
5^ ', are three basic channels for lepton scattering on nuclei, namely elastic, quasielastic, and inelastic 
processes. In the case of elastic scattering (z/ = Q'^/2Ma, where Ma is nuclear mass) the leptons 
are scattered off the nucleus leaving the nucleus in its ground state. Quasielastic scattering 
(z/ ~ Q^/2M, where M is nucleon mass) corresponds roughly to direct collisions with the 
individual nucleons inside the nuclei. The inelastic scattering occurs when the pion threshold 
is reached (z/ > Q'^/2M + m^ where m.,^ is the pion mass). At the Born level both Q^ and 
1/ are fixed completely by measuring the scattering angle and the energy (momentum) of the 
scattered lepton. However, at the level of radiative corrections, in case of presence of real 
radiated photons, the fixation is removed and the four-momentum of the radiated photon has 
to be included in the kinematical variable calculation. Such elastic, quasielastic and inelastic 
processes with radiation of a real photon are known as radiative tails from the elastic {(Tei) 
and the quasielastic ((Xg) peaks and from the continuous spectrum (crj„) called hereafter shortly 
elastic, quasielastic and inelastic radiative tail. 

The total radiative correction at the lowest order is obtained as the sum of these contri- 
butions together with loop corrections (o"„) coming from effects of vacuum polarization and 
exchange of an additional virtual photon: 








a) b) c) d) e) 

Figure 1: Feynman diagrams contributing to the Born and the radiative correction cross 
sections in lepton-nucleus scattering. 
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In this report we present a Monte-Carlo generator for events with a possible radiation of 
a real photon in deep inelastic lepton scattering (DIS) on polarized and unpolarized targets. 
It should be noted that only the single photon exchange (no Z^ exchange) and only pure 
QED corrections are considered. Hence the use of the generator is at present restricted to 
experiments where the electroweak contributions and corrections are small, as e.g. in the 
experiment HERMES at DESY. The events are generated in accordance with their contribution 
to the observed total cross section given by 

(Tabs = (Tnon-radi^) + (^ini^) + (Tg + ael- (2) 

The first term (7non-rad{^) contains not only the contribution from the Born process but also 
the contributions from loop corrections (a^) and from multiple soft photon production with a 
total energy not exceeding a cut-off parameter A. In handling the radiative corrections two 
approaches have been used, the one developed by Mo and Tsai [|l| and the one given by Bardin 
and Shumeiko [0]. These two formalisms and the significance of the cut-off parameter A will be 
discussed below in more detail. The generator contains two patches which can be considered 
as independent generators. The first one is based on the FORTRAN code POLRAD 2.0 0] 
handling the polarized lepton- nucleon scattering. The other is based on FERRAD [Q, a code 
that deals with the unpolarized case. They will be referred to as the POLRAD and FERRAD 
generator throughout the text. A detailed analytical and numerical comparison of the results 
given by the two generators is performed. 

2 The Kinematics and Main Stages of Generation 

The Feynman diagram for the Born (or one photon exchange) process is presented in Fig.|I]a. 
Its kinematics is completely defined by the scattering angle 6 and the energy E' of the scattered 
lepton, two variables which are usually measured. All other inclusive kinematical variables can 
be expressed in terms of 6 and E'. Having in mind the laboratory frame, i.e. the frame with 
the target nucleus at rest the relevant variables are 



Q^ = AEE' sin^ -, u = E-E', x 



2' ' 2MZ/' .gS 



where E is the beam energy, Q^ is the negative of the four-niomentum transfer squared, z/ the 
energy transfer, x the Bjorken scahng variable, y the normahzed energy transfer, and W"^ the 
mass squared of the hadronic final state. Radiative QED corrections (RC) at the lowest order 
are described by the set of the Feynman graphs shown in fig. P^-e. Contributions to the lowest 
order RC cross section come from the diagrams b) and c) and from the interference of the loop 
diagrams d) and e) with the Born matrix element a). 

An event registered in the detector with certain values of E' and 6 for the scattered lepton 
can either be a non-radiative or a radiative event, i.e. an event containing a real hard radiated 
photon. For a radiative event there are apart from E' and 6 three additional variables necessary 
to fix the kinematics of the radiated real photon. A possible choice is the photon energy E^ 
and the two angles 6^ and 0^, where 6^ is the angle between the real and the virtual photon 
momenta k and q = ki — k2 {ki^2 being the momenta of the initial and the scattered lepton) and 
(f).y is the angle between the planes defined by the momenta {ki, ^2) and {k, q). For events with 
the radiation of a real photon the kinematical variables describing the virtual photon and used 
to generate the hadronic final state differ from eq.(|^) because the substitution q ^ q — k has 
to be made in their definition. The variables obtained after this substitution will be referred 
to as 'true' ones: 



W, 



true 



W^ - 2E^{iy + M - Ju^ + Q^ cos^^), Utme = v - E. 
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QLe = Q' + ^E,{iy - Jz/2 + Q^cose. 



'1)1 



Xij- 



Q 



(4) 



true 



2Mu 



true 



For non-radiative events the true kinematics exactly coincide with (|^). 

In case of the elastic and quasielastic radiative processes the mass squared of the hadronic 
final state is fixed imposing additional constraints on the true kinematical variables. Indeed 
the following ranges are allowed: 



X < Xtrue < 1 for (Ti„ 



Xtrue = 1 for (To 



Xtrue = Ma/M for a el 



^min — ^true — W 
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(5) 



where 



and 
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,2{l-Xr){l±VTT^)+l' 



7^ + 4Xj.(l — Xr 



Y = 4:M^XyQ\ Xr = X/Xtr 



(6) 
(7) 



The Monte-Carlo procedure for the generation of events with a possible photon radiation is 
as follows: 

It starts off with the generation of the kinematics of the scattered lepton and the calculation of 
an event weight from these kinematics. Then the appropriate scattering channel (non-radiative; 
elastic, quasielastic or inelastic radiative tail) has to be chosen according to their contribution 
to the total observed cross section (see (0)). If a radiative channel is selected the radiated 
photon has to be generated and the values of the kinematical variables have to be re-computed 



in order to obtain the true values. For each event the weight has to be recalculated. The new 
weight is defined as the ratio of the radiatively corrected and the Born cross section. After this 
recalculation the weighted sum of all events (generated originally in accordance with the Born 
cross section) gives the observed cross section. 

The recalculation of the weight requires the knowledge of the cross section integrated over 
the photon momentum. This is done differently in the two theoretical approaches which will be 
discussed in detail in the subsections below. In the Mo-Tsai approach an additional parameter 
A is introduced dividing the integration region into a soft and a hard photonic part. No such 
parameter is used in the Bardin-Shumeiko approach. However, in both approaches a minimal 
photon energy -Emm is adapted in generating a radiated photon. Only if the photon energy is 
above this value the kinematics of the photon is calculated and the event becomes a radiative 
one. The actual value of -Emm depends on the aim of the photon generator. In general photons 
should be detected in the calorimeter. So the energy threshold of the calorimeter sets the value 
for E'^ ■ 

In addition one needs spin averaged and spin dependent structure functions, quasielastic 
suppression factors and elastic formfactors to calculate the Born and observed cross sections. 
This is taken either from fits to experimental data or from model predictions. A recent 8- 
parameter fit from ref.[B] for the proton and deuterium structure function -Ff' [x, Q'^) and their 
combination for ^He (-^2^'' = 3-^2'+ 3-^2^) has been used. The spin dependent structure function 
gi{x, Q"^) is constructed from a fit to the spin asymmetry Ai [Q according to 

g,{x,Q')=Mx)F^{x,Q') = ^^iMZ^M!) (1 + ^2)_ (g) 



2x(l + R{x,Q^' 



For the structure function R{x, Q"^) being the ratio of longitudinal to transverse virtual 
photon absorption cross sections Whitlow's fit was used and assumed to be A-independent. 
No contributions are assumed from the spin dependent structure function g2{x,Q'^)- For the 
elastic formfactors a fit to experimental data [|^ is utilized. The quasielastic radiative tail is 
calculated in accordance with ref. @]. In contrast to the elastic case where the treatment is exact 
(within the considered order of perturbative theory) here the Y-scaling hypothesis and the 
peaking approximation (quasielastic structure functions are evaluated at the quasielastic peak) 
are applied. 

The two Monte Carlo generators for radiative events available in RADGEN for the unpolar- 
ized and polarized cases are described below. The theoretical approaches, the approximations 
made, and the dependence of the radiative corrections on the artificial parameter A and on the 
value of the photon detector threshold E^^in are discussed. 

2.1 The generator for the unpolarized case 

The unpolarized generator is based on the FORTRAN code FERRAD35 . The code calculates 
the radiative correction to deep inelastic scattering of unpolarized particles in accordance with 
the analytical formulae given by Mo and Tsai [|l|. 

The result for the lowest order radiative correction of the cross section is 

a = 5r{A){1 + 5yert + Svac + ^sm)^"!^, + ffe/ + CTg + ai„(A), (9) 
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where o"i^ is the one photon exchange Born cross section and S^aci ^vert, and Sgm are corrections 
due to vacuum polarization by electron and muon pairs, vertex corrections and residuum of the 
cancellation of infrared divergent terms independent of A: 
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where m is the lepton mass and Li2(x) 
function. 



(10) 



Jq dyln{l — y)/y is the dilogarithm or Spence 



The cross sections aei, Cq, and ain are the contributions from radiative processes (i.e. in- 
cluding real photon radiation) for elastic, quasielastic and deep inelastic scattering. They can 
be calculated in terms of the radiative tail from the jt/i mass level aj 



(Tei = Crez(M4,l), (Tg = (Xg^M, 1) , 



W 

a.n= J dM,aUM,,e. 



(11) 



where Mh = JlVtme and aj {j = el, q, in) has the form of an integral over 6^, the angle between 
the real and the virtual photons: 



"max 

a,{M,ema.)= j d9^ ToiWiiT, + T2 + Ts + T^ + T,) + W(in + Tr)). 



(12) 



The explicit form of the structure functions W12 depends on the type of the tail. The 
terms Tj are kinematical factors (see ref. |jl| for details). The integration limit 6max is defined 
as 6max = fnin(7T, 6/^) where ^a can be found from the relation 



W, 



true 



W - 2A(z/ + M-Jiy^ + Q^ cos(0a)) 



(13) 



An artificial parameter A had to be introduced to divide the integration region over the 
photon energy into two parts, the soft and the hard energy regions. The hard energy region 
(Tj„(A) can be calculated without any approximations. The soft photon part is calculated 
for photon energies approaching zero. After cancellation of the infrared divergences and a 
resummation of soft multiphoton effects the correction factor is given by 



Sr{A) = exp 



a / E E' 

- In— + ln — 
7^ \ A A 




(14) 



Note that we follow ref.|l| in fixing that part of the soft photon correction which should 
be exponentiated. Another possible variant is the well-known prescription of ref.[0. The 
dependence on A is the deficiency of the method. On the one hand this artificial parameter has 
to be chosen sufficiently large to avoid numerical instabilities because for A = the contribution 
o"j„(A) gets infinite. On the other hand it has to be chosen sufficiently small to reduce the soft 
photon region which is calculated only approximately. If the threshold over the photon energy 
is low enough, it is customary to use the same value for the parameters A and El 
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2.2 The generator for the polarized case 

The polarized generator is constructed utilizing the FORTRAN code POLRAD 2.0 [^, |TT|. 
This code is based on the method of covariant cancellation of infrared divergences developed 
by Bardin and Shumeiko |0. The final formula for the cross section free of infrared divergences 
does not include any artificial parameter like A. However, -E^jn has to be introduced if the 
code is used considering the generation of real radiated photons. The cross section formula of 



ref.|TT| can be rewritten into a form similar to eq.( 



Cr = 6R{E^^i,^){l + 6^ert + ^vac + ^sm)(^l^ + (^add{El,in) + (^el + ^-g + Crm(-Emm)- (15) 

In the ultrarelativistic approximation (Q^ ^ m?) the corrections 5yerti ^vaa and 5sm. take 
the same form as in eq.(|10|) Q. The covariant form of the radiative tails for a nuclear target 
with atomic number A is as follows: 






' mm 
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3 TAya-^ 4 1^^ 9A/f2pi-2 

The integration limits are determined by Tmax,min = Q^/2M^x(l ± yl + 4M2x2/(5^), Rmin = 
2ME'^i^, and -Rmax = iW^ — (M + m7r)^)/(l + r). The quantities % are kinematical factors, 
and 3j are some combinations of the DIS structure functions, elastic formfactors or quasielastic 
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response functions (see ref. |jTl| for details). The sum over i corresponds to contributions from 
spin averaged {i = 1,2) and spin dependent {i = 3,4) structure functions. In the case of the 
elastic and quasielastic radiative tail there is no integration over R because the R range is 
reduced to the constant values R = Rei a and R = R^i {Rei = (2Mz/ — Q'^)/{1 + r)) due to 
the fixed final hadron mass. Invariants with the index "A" contain the nucleus momentum p^ 
instead oi p {p\ = M|). 

Apart from the usual ultrarelativistic approximation another one was made when the eqs. 
d^) and (|1^) were obtained, namely the photon energy was considered to be small E^ -C E, E' 
in the region E^ < -E^j„. In POLRAD, the term aadd{E'^^^) is added to take the difference of 
this approximation to the exact formula into account: 

2ME'' - 

dR 



(Jadd{E^in) = -tt^y j dT^hiiir) J 



i=l I 



k. 



R 




^?/«M / ^^w^TW*-'"''"' ■ '"' 



Therefore, the results obtained with the help of POLRAD are independent of the threshold 
parameter -E^m by construction. 

It should be noted that contributions from r-leptons and quark loops were included into 5yac P^ - 
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3 Comparison and Numerical Results 

3.1 Analytical comparison 

The expressions for the spin independent part of the radiative tails in eq. (0) have been 
compared analytically with the radiative tails given by eqs.(|TT]) and (|T^ for the unpolarized 
case using the algebraic programming system REDUCE 3.5. It turned out that 



3 



-2M^v^ + Q^Y.^u{t)R^-^ = n + Tr. (18) 

i=i 

From these relations one can easily show that the formulae for the radiative tails are exactly 
the same when no spin is considered. 

All corresponding terms in (Rf) and (p!5D coincide exactly and there is only one additional 



term aadd in eq. (p!5D . This term can be considered as characterizing the quality of the soft 
photon approximation and can be calculated analytically by expanding it in the parameter 

e = 2MEIJQ': 



aaddiElJ = e-ao - e^{SrF^ + S22F2) + 0{e'), (19) 



where ctq is obtained in terms of the Born cross section cxo: 



,„,^^_,J_2F. + .^ + Q=^^- 







dx ' "^ ^ ) "' \ ^ ' "^ dx ' ^ dQ'^ 
The quantities Jq, 6r, and 622 depend only on kinematical variables: 

I — y \ m J \ m J m^ 



(21) 



522 = -4^-^ ln(l -y)- 2 Jo - ^5.- 

y Q 

The first term in eq.(|19D corresponds to the approximation ii^^ = in arguments of structure 
functions (terms with derivatives) and the photonic propagator (terms without derivative). The 
second term in eq.(|l9|) appears when the approximation is applied to kinematical terms like 
T1..7 as in eq . ([T^) . 

The numerical analysis performed for HERMES kinematics shows that the difference be- 
tween both generators in the unpolarized case due to this additional contribution is negligible. 
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Figure 2: The weighted Xtme o-nd Qfj-ue distributions as generated by FERRAD (histogram) 
and POLRAD with polarization part switched off (points) a), b); and by POLRAD with parallel 
(histogram) and antiparallel (points) spin configurations c), d). 



3.2 Numerical comparison 

In order to compare these two generators numerically Monte-Carlo event samples of 100k events 
each have been generated for a '^He target and kinematical cuts relevant for the HERMES ex- 
periment 11 have been applied (Q^ >i GeV^, W^ >4 GeV^, y <0.85, 0.037< 6 <0.14rad). 
Weighted distributions (A^^) of Xtr-ue and Ql^-ue generated by the FERRAD and the POLRAD 
generators are shown in FigJ^. In Figs. ^ and ^b the unpolarized case is compared, i.e. FER- 
RAD and POLRAD with the polarization part switched off. As can be seen the distributions 
are in agreement. It should be noted that the events at x ~ 1 and x ~ 3 correspond to the 
quasielastic and elastic radiative tails, respectively. In Fig. |^) and ^) a comparison of the 
parallel and antiparallel spin configuration in POLRAD is given. Here a systematic difference 
is seen which is expected from a spin asymmetry different from zero at small values of x. 

The distributions of the radiation angles 6^ and 0^ as defined in the so-called Tsai-system 
in which the z-axis is along the direction of the virtual photon and the y-axis is normal to the 
scattering plane (see ref. [Q]) and of the energy E^ are displayed in Fig.||. These distributions 
are generated for a certain point in the kinematical plane of the scattered lepton, i.e. for x 
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Figure 3: The distribution of the radiation angles 9^ a) and 0^ h) and of the energy c) of the 
radiated photon for x = 0.1 and y = 0.8. The two-dimensional distribution d) shows 9^ vs 0^. 



= 0.1 and y = 0.8. The two peaks seen in the 9^ distribution of Fig.^ correspond to the s- 
and p-peak emerging from colhnear radiation along the direction of the incoming and outgoing 
lepton, respectively. 

Further results on the numerical comparison of the polarized and unpolarized case in the 
covariant and Mo-Tsai approach, respectively, can be found in the refs. [|^, |13]. Fig.4 of ref. [|13 



presents the A-dependence of the radiative correction factor obtained by the code FERRAD. 
The general conclusion was that the best choice is A equal to about 0.1% of the beam energy. 
In this case the results should not depend on A. 

The other task of the generators is to generate a radiative correction factor (Tobs/o'Born 
necessary to recalculate the weight of each event. A typical distribution of the RC factor is 
shown in Fig.^. Apart from this integration characteristics it is useful to consider also the 
behavior of the RC factor as a function of kinematical variables. Both the observed and 
Born cross sections can be split into spin averaged and spin dependent parts (^(j°bs,Born _ 
^obs, Born _^p^p^^obs, Born ^ Pb,t are beam and target polarizations) allowing to define the RC factors 
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Figure 4: The typical distribution of 
the radiative correction factor. 
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Figure 5: The spin averaged (full 
lines) and spin dependent (dashed 
lines) RC factor as a function of y for 
for different values of x. 



for the unpolarized and polarized case separately: 
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(22) 



Here we keep in mind that there is a region (x ~0.4 for ^He case) where cr^°'''"'=0. The y- 
dependence of these factors for several values of x is plotted in Fig.|^. We would like to mention 
that the RC to the polarization asymmetry Aobs,Born = cr°*^'^'"""/cr°*'*'^°'"" can be written in 
terms of 6u,p as 

A A XX 

(23) 
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As can be seen from Fig.|^ the correction factors 6p and 6u have similar behaviour and become 
almost equal with increasing x. The large corrections coming from the factorized part of eqs.(^ 
and (^) (the first term) cancel exactly in the numerator of (p3[) and only the unfactorized part 
of the corrections coming from the radiative tails contribute to the RC of the polarization 
asymmetry. This is the reason why there is a large correction to the total cross section and 
only a small one to the asymmetry. 



4 Conclusion 

The Monte-Carlo generator RADGEN 1.0 including radiative events in DIS on polarized and 
unpolarized targets have been presented. Two patches of the generator, one for the unpolarized 
case based on the FORTRAN code FERRAD35 and the other for the polarized case based on 
the FORTRAN code POLRAD 2.0 were compared and tested. This generator is implemented 
into the HERMES Monte-Carlo program HMC. The FORTRAN code of the Monte Carlo 
generator RADGEN 1.0 is available on request. 
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Appendix: Main options and keys 

In this Appendix we briefly discuss options and keys which should be specified by a user in 
order to run RADGEN. They all are gathered in one input file 'input.dat' (default values are 
given in brackets): 

• nevent [1000] - number of events; 

• ixytest [0] - usage of look-up table for channel generation. The information about results 
of RC calculation in the 50x50 points over x and y is stored in this look-up table. The 
key defines the way how to use this information: 

1 two-dimensional spline 

2 PINT interpolation 

generation of the look-up table 
-1 No use of the look-up table 

• ige [2] - unpolarized (1) or polarized (2) generator is used; 

• ire [1] - type of target: proton ( 1), deuteron (2), hehum-3 (3); 

• el [27.5 GeV] - beam energy; 

• pi and pn [0 and 0] - degrees of beam and target polarization; 

• ikey [0] - type of generation of leptonic kinematics 

1 — X and y are generated in accordance of the cross section within kinematical cuts, 
— X and y are read from the standard input; 

• ntk [35] - number of bins in each of the seven parts of the integrand over the azimuthal 
photonic angle. The region over the azimuthal photonic angle is divided into seven parts 
in order to locate the peaks. This key defines the number of bins in each part of the 
integrand over the azimuthal photonic angle where the cross section is calculated. 

• nrr [200] - number of bins in the integrand over photonic energy. This key defines the 
number of the bins over photon energy where cross section of the process is calculated 
and stored. 

• demin [0.1 GeV] - A the threshold of photonic energy. Events with photonic energy 
smaller than A are considered as non-radiative. 

• Kinematical cuts on y, x, Q^, W^, 9. 
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